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Abstract 

This paper describes the first measurement of ^-quark fragmentation fractions into bottom 
hadrons in Run II of the Tevatron CoIUder at Fermilab. The result is based on a 360 pb^^ sam- 
ple of data collected with the CDF II detector in pp collisions at ^/s = 1.96 TeV. Semileptonic 
decays of B^, B~ , and mesons, as well as A|J baryons, are reconstructed. For an effective 
bottom hadron pT threshold of 7 GeV/c, the fragmentation fractions are measured to be fu/fd = 
1.054 ± 0.018 (stat) t'^ofi (sys) ±0.058(6), /,/(/„ + /d) = 0.160 ±0.005 (stat) tO:^^! (sys) (g), 
and fAjifa + fd) = 0.281 ± 0.012 (stat) t^^^^ (sys) t^^^^ {B), where the uncertainty B is due to 
uncertainties on measured branching ratios. The value of fsjifu ± fd) agrees within one standard 
deviation with previous CDF measurements and the world average of this quantity, which is dom- 
inated by LEP measurements. However, the ratio fh^^/{fu ± fd) is approximately twice the value 
previously measured at LEP. The approximately 2 a discrepancy is examined in terms of kinematic 
differences between the two production environments. 

PACS numbers: 13.20.He, 13.30.Ce, 14.20.Mr, 14.40.Nd, 14.65.Fy 
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I. INTRODUCTION 



Bottom quarks, b, produced in pp collisions combine with anti-quarks or di-quarks to 
form bottom hadrons. In this process, called fragmentation, the color force field creates 
quark- ant iquark pairs qq that combine with the bottom quark to create a B meson \bq) or 
b baryon \bqiq2)- Since the fragmentation process, which is governed by the strong force, 
cannot be reliably calculated by perturbative QCD [1-3], the fragmentation properties of 
b quarks must be determined empirically. This paper describes a measurement of the species 
dependence of the 6-quark fragmentation rates into bottom hadrons produced in pp collisions 
at center of mass energy ^/s — 1.96 TeV during Run II of the Tevatron collider at Fermilab. 

The probabilities that the fragmentation of a 6 quark will result in a B~ \bu), \bd), or 
Bg \bs) meson or a A° \bdu) baryon are denoted by /„, fd, fs, and /a^, respectively. In this 
paper, fg indicates the fragmentation fraction integrated above the momentum threshold of 
sensitivity in the data: fq = fq{pT{B) > pjp'^'") [4]. In the case that the fragmentation frac- 
tions are momentum dependent, the measured fragmentation fractions are proportional to 
the relative yields of the bottom hadrons integrated above the effective p™". The contribu- 
tions from the production of excited bottom hadrons that decay into final states containing 
a B~ , Bg meson or baryon are implicitly included in this definition of the fragmenta- 
tion fractions, fq = B{b BqX). Throughout the paper, unless otherwise noted, references 
to a specific charge state arc meant to imply the charge conjugate state as well. 

In Run I of the Fermilab Tevatron, which collected data from 1992-1996, the fraction 
of Bg mesons produced relative to the number of 5° mesons was measured 2 cr higher 
at CDF [5-7] than at the LEP experiments [8-10]. Interestingly, the time-integrated flavor 
averaged mixing parameter, x — fdXd + fsXs, where Xd and Xs are the time- integrated 
mixing parameters of B^ and B^ mesons respectively, was also measured 2 cr higher in 
Run I [11, 12] than the LEP averages of the same quantity [13-18]. This second discrepancy 
led to speculations about possible sources of the enhanced average mixing rate at a hadron 
collider relative to electron-positron collisions, including suggestions that new physics may 
be the source of the disagreement [19]. Since the average momentum of b quarks produced at 
LEP, {p{b) ) ~40 GeV/c, is significantly higher than at the Tevatron, {p{b) ) ~ 10 GeV/c, it 
is also possible that the fragmentation process depends on the 6-quark momentum. Another 
possible explanation is that fs is higher at the Tevatron than at LEP due to the different 
initial mechanism of 6-quark production. Of course, a more mundane possibility is that the 
Run I results relating to fs are simply statistical ffuctuations. To shed light on the question of 
whether 6-quark fragmentation is different in a hadron environment than in e'^e" collisions, 
the fragmentation fractions are measured in CDF Run II with high statistical precision and 
an updated treatment of the lepton-charm sample composition. 

The analysis strategy is as follows. Semileptonic decays of bottom hadrons, B — > 
i~'p£DX, where £~ stands for electron or muon, and D represents a charm meson or baryon, 
in case of semileptonic bottom baryon decays, unless otherwise specified, provide large sam- 
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pies for studying the fragmentation properties of b quarks. This measurement determines 
the 6-quark fragmentation fractions by reconstructing five semileptonic signatures, £~D'^, 
£~D^, £~D*~^, £~Df, and ^~A+. The selection requirements are kept similar among the 
five lepton-charm channels in order to cancel as many systematic uncertainties as possible. 
The final signal requirements, though similar, have been selected to maintain good accep- 
tance for the individual decays, which have different kinematic features. The reconstructed 
£^D signal yields, originating from the various B £^UiDX semileptonic decays, arc then 
related to the numbers of bottom hadrons {B~ , 5°, B^, or A^) produced in the 6-quark frag- 
mentation process. Since the neutrino from the semileptonic bottom hadron decay is not 
reconstructed, the missing energy in the decay allows semileptonic bottom hadron decays 
to excited charm states to contribute to the five final state decay signatures. This results in 
"cross-talk" between the bottom hadron species, particularly between the B mesons. The 
observed semileptonic £"0 decay signatures are related to their corresponding parent bot- 
tom hadrons through a procedure used to extract the sample composition, as described later 
in the text. In order to reduce systematic uncertainties in trigger and tracking efficiencies, 
the 6-quark fragmentation fractions are measured relative to fd- This means that the rel- 
ative fragmentation fractions fu/ fd, fs/{fu + fd) and /aj/(/u + fd) are extracted from the 
five lepton-charm yields, taking the sample composition into account. Since the fragmen- 
tation of b quarks into b baryons other than the A|| are ignored, a constraint requiring the 
fragmentation fractions /„, f^, fs, and f\^ to sum to unity is not applied. 

This paper is organized as follows. The semileptonic signal reconstruction is discussed 
in Section 11, while the sample composition procedure used to relate the lepton-charm sig- 
natures to the parent bottom hadron is described in Section III. The efficiencies needed 
to extract the sample composition are determined in Section IV. The fit to the fragmenta- 
tion fractions is detailed in Section V. Finally, the systematic uncertainties assigned to the 
measurement are described in Section VI and the final results are discussed in Section VII. 

II. DATA RECONSTRUCTION 
A. Experimental Apparatus 

The data used in this measurement represent an integrated luminosity of approximately 
360 pb^^ collected with the CDF II detector between February 2002, and August 2004. 
The CDF detector employs a cylindrical geometry around the pp interaction region with the 
proton direction defining the positive 2;-direction. Most of the quantities used for candidate 
selection are measured in the plane transverse to the 2;-axis. In the CDF coordinate system, 
(f is the azimuthal angle, 9 is the polar angle measured from the proton direction, and 
r is the radius perpendicular to the beam axis. The pscudorapidity rj is defined as ?7 = 
— ln[tan(^/2) ]. The transverse momentum, pt, is the component of the track momentum, p, 
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transverse to the z-axis {pr — p-sin 9), while Et — E-sin 9, with E being the energy measured 
in the calorimeter. 

The CDF 11 detector features excellent Icpton identification and charged particle tracking 
and is described in detail elsewhere [20, 21]. The parts of the detector relevant to the 
reconstruction of scmileptonic bottom hadron decays used in this measurement are briefly 
summarized below. The detector nearest to the pp interaction region is a silicon vertex 
detector (SVXII) [22], which consists of five concentric layers of double-sided sensors located 
at radii between 2.5 and 10.6 cm. An additional single layer of silicon (LOO) [23] is mounted 
on the beam pipe at radius r ~ 1.5 cm, but the information from this detector is not used 
in this measurement. In addition, two forward layers plus one central layer of double sided 
siheon located outside the SVX at radii of 20-29 cm make up the intermediate sihcon layers 
(ISL) [24]. Together with the SVXII, the ISL detector extends the sensitive region of 
the CDF II tracking detector to I77I < 2.0. CDF's sihcon system provides three-dimensional 
track reconstruction and is used to identify displaced vertices associated with bottom hadron 
decays. The measurement of the momentum of charged particles in the silicon detector is 
significantly improved with the central outer tracker (COT) [25], an open-cell drift chamber 
with 30,200 sense wires arranged in 96 layers combined into foiir axial and four stereo super- 
layers (SL). It provides tracking from a radius of ~40 cm out to a radius of 132 cm covering 
\z\ < 155 cm. The track reconstruction efficiency of the COT is found to be {99.Q^Q'g)% 
for charged particles with Pt > 1.5 GeV/c [26] and > 94% [27] for charged particles with 
Pt = 0.4 GcV/c. For high-momentum charged particles, the Pt resolution is found to 
be (j{pt)/pt = 0.0015 pr/GcV/c. The COT also provides specific energy loss, dE/dx, 
information for charged particle identification with a separation between pions and kaons 
of approximately 1.4 cr [28]. The central tracking system is immersed in a superconducting 
solenoid that provides a 1.4 T axial magnetic field. 

Electromagnetic (CEM) [29] and hadronic (CHA) [30] calorimeters are located out- 
side the COT and the solenoid, where they are arranged in a projective-tower geometry. 
The electromagnetic and hadronic calorimeters are lead-scintillator and iron-scintillator 
samphng devices, respectively. The energy resolution for the CDF central calorimeter 
is (t{Et)/Et = [(13.5%//&;)^ + (1.5%)2]V2 for electromagnetic showers [29, 31] and 
(t{Et)/Et = [{75%/^/E^f + (3%)2]V^ for hadrons [21, 30], where Et is measured in GeV. 
A layer of proportional chambers (CES), with wire and strip readout, is located six radiation 
lengths deep in the CEM calorimeters, near the electromagnetic shower maximum. The CES 
provides a measurement of electromagnetic shower profiles in both the (/?- and z-directions 
for use in electron identification. Muon candidates are identified with two sets of multi-layer 
drift chambers and scintillator counters [32, 33], one located outside the calorimeters (CMU) 
and the other (CMP) behind an additional 60 cm of iron shielding, equivalent to approxi- 
mately 3 pion interaction lengths. The CMU provides coverage for particles with \r]\ < 0.6 
and Pt > 1-4 GeV/c. The CMP covers the same pseudorapidity region, but identifies muons 
with Pt > 2.0 GeV/c with higher purity than muons reconstructed in the CMU only. 
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B. Trigger Requirements 

CDF uses a three-level trigger system [21], where each level provides a rate reduction 
sufficient to allow for processing at the next level with minimal dead-time. At level 1, data 
from every beam crossing are stored in a pipeline memory capable of buffering data for 
~ 5.5/i.s. The level 1 trigger either rejects an event or copies the data into one of four level 2 
buffers. At level 2, a substantial fraction of the event data is available for analysis by the 
dedicated trigger processors. Events that pass the level 1 and level 2 trigger selection criteria 
are then sent to the level 3 trigger [34, 35], a cluster of computers running a speed-optimized 
reconstruction code. Events selected by level 3 are written to permanent mass storage. 

Tracking plays a significant role in the triggers utilized for this analysis. Semileptonic 
B — > £~P^DX decays are recorded using a trigger that requires a lepton and a track displaced 
from the interaction point and identified with the sihcon vertex trigger (SVT) [36]. The 
decay topology of semileptonic B decays is sketched in Fig. 1. Tracks are reconstructed at 
level 1 with the extremely fast tracker (XFT) [37] by examining COT hits from the four 
axial super-layers. The XFT provides r-ip tracking information and can identify tracks with 
Pt > 1-5 GeV/c with high efficiency (> 90%) and good transverse momentum resolution, 
cr{px)/pT = 0.016 pT/[GeV /c]. XFT tracks can be matched with cither calorimeter clusters 
to identify electron candidates or with track segments in the muon detectors to identify 
muon candidates. The XFT tracks are extrapolated into the silicon detector system, where 
the SVT uses the SVXII measurements of charge deposits from charged particles to form 
simplified tracks. In addition, the SVT determines the distance of closest approach in the 
transverse plane, do, with respect to the pp beam line, which is determined from a time- 
dependent line fit to the locus of primary interaction vertices determined from all tracks 
available at trigger level (see Fig. 1). The impact parameter resolution of the SVT is 
approximately 50 /im [36, 38], which includes a contribution of 35 /im from the width of the 
pp interaction region [39]. 

The primary trigger used in this measurement requires that the lepton and the dis- 
placed track (SVT track) must have transverse momentum values greater than 4 GcV/c and 
2 GcV/c, respectively. The displaced track's impact parameter, do, must exceed 120 //m 
and be less than 1 mm to reject decay products of long-lived hadrons decays such as 
Kg or A°. The opening angle, A0, between the lepton and SVT track is required to 
satisfy 2° < A0(£^, SVT track) < 90° to increase the probability that the two tracks 
originate from the same B hadron. Additionally, the invariant mass between the trig- 
ger lepton and SVT triggered track must be less than the nominal bottom hadron mass, 
m{i-, SVT track) = ^(pM(£-) + p'^(SVT))2 < 5 GeV/c^ where the SVT track is assumed 
to have the pion mass. The trigger lepton requirements are described in conjunction with 
their analysis selections in Section II CI. Events that pass these trigger requirements are 
recorded to the lepton plus SVT trigger data stream for further analysis. In this measure- 
ment both the muon and electron plus SVT trigger data (e-|-SVT and /x-l-SVT) are used. 
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FIG. 1: Sketch of semileptonic -B-decay topology in the transverse plane, where Vp is the primary 
vertex, Vb is the decay vertex of the bottom hadron, Vd is the decay vertex of the charm hadron, 
and do is defined in the text. "SVT" indicates the track selected by the displaced track SVT 
trigger, which is also defined in the text. 

An additional trigger utilized for selecting B events is the two-track trigger (TTT), which 
requires two displaced tracks. Large semileptonic 13 samples are also available with this 
trigger [40, 41], although the false lepton background is much larger as well. Semileptonic 
events from the TTT are used in this analysis for a study of the systematic uncertainty 
arising from false leptons. 

C. Data Selection and Reconstruction 

Events from the lepton plus SVT trigger data stream are used to reconstruct semileptonic 
bottom hadron decays in this analysis. First, trigger leptons are identified by re-confirming 
the trigger decision with offline quantities after event reconstruction. Charm candidates are 
then reconstructed, with the SVT track required to match one of the daughter tracks from the 
charm decay. The selections on the lepton-charm signals obtained are optimized to reduce 
combinatoric background and improve signal significance. Non-combinatoric backgrounds 
in the charm signals are handled separately. 

1. Trigger Lepton Identification 

The data analysis begins by identifying the trigger leptons from the e-|-SVT and /i+SVT 
trigger streams. The electron candidates are identified by requiring the following selection 
criteria. The longitudinal shower profile must be consistent with that of an electron shower, 
with a leakage energy from the CEM into the CHA of less than 12.5%, in order to sup- 
press hadron contamination. The lateral shower profile of the CEM cluster is required to 
be consistent with a profile obtained from test beam electrons after appropriate corrections. 
The association of a single track with the calorimeter shower is made based on the position 
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matching at the CES plane, with both IA2; • sin^| < 5 cm and rlA^?! < 3 cm conditions 
required. To achieve good agreement between data and Monte Carlo (MC) simulation (see 
Sec. IV A), an isolation requirement is applied to the trigger electron candidates by requir- 
ing that exactly only one track is found that projects to the CEM towers used to define 
the electron energy. To reconfirm electron trigger cuts, the offline reconstructed Et and 
Pt of the electron candidate are required to be greater than 4 GeV and 4 GeV/c, respec- 
tively. Additionally, electron candidates from photon conversions in the detector material 
are removed by rejecting those electron candidates that have a small opening angle with 
oppositely charged particles in the event. 

Trigger muon candidates are reconstructed by extrapolating tracks measured in the COT 
to the muon system, where they are matched to track segments (stubs) reconstructed in the 
muon chambers. A CMU or CMP stub is required to have hits in at least three out of the 
four layers of planar drift chambers. Trigger muons are required to have hits in both the 
CMU and CMP muon chambers. The separation between a track segment reconstructed in 
the muon chamber and the extrapolated COT track is computed. The uncertainty in this 
quantity is dominated by multiple scattering in the traversed detector material. For good 
track to stub matching, this separation is required to be less than 15 cm and 20 cm in the 
r(p-yiew for CMU and CMP, respectively. The transverse momentum of a muon candidate 
reconstructed offline is required to be greater than 4 GeV/c. 

2. Charm Candidate Selection 

The SVT track is required to match one of the final state tracks in the five reconstructed 
charm signals: D° ^ K-n+, D*+ (^ K~n+)n+, D+ K-n+n+, D+ ^ (^ 

K~^K~) 7r~^, and A^ pK~Ti^. Only well-reconstructed tracks with pr > 0.4 GeV/c and 
at least three silicon r-(p hits are retained for offline analysis. To ensure good track quality, 
all charm daughter tracks, except for the soft pion from the D*^ decay, are required to 
have at least five hits in at least two axial and two stereo COT super-layers. There are no 
COT requirements on the D*+ soft pion. During data reconstruction the track parameters 
are corrected for the ionization energy loss appropriate to the mass hypothesis under con- 
sideration. In addition, tracks are required to be fiducial in the COT, so that only tracks 
which are well-described by the simulation (see Sec. IV A) are used for further analysis. In 
particular, tracks that fall within \z\ < 1.5 cm of the COT mid-plane, where no track infor- 
mation is recorded, and tracks that originate outside of the COT volume at \z\ > 155 cm 
are excluded from the analysis. In addition, all tracks must at least pass through the axial 
SL 6 before exiting the COT. This means the exit radius of the track must be greater than 
the radius of the sixth super-layer rsLg = 106 cm. This requirement is tightened for the 
SVT trigger track and the trigger Icpton. Both tracks must pass through SL 8 of the COT 
{t^SLs = 131 cm) as required in the trigger. The invariant mass of the — > K~7t~^ and 



12 



D+ K-t:^t:+ is reconstructed within [1.40,2.00] GeV/c^and [1.70,2.00] GeV/c^, respec- 
tively. The reconstructed — > ^vr"'" mass is required to be within [1.75, 2.2] GeV/c^, while 
the A+ pK~7r^ is reconstructed within [2.15,2.40] GeV/c^. Finally, the reconstructed 
charm signals are combined with the triggered lepton in a three-dimensional kinematic fit 
constraining all tracks to a common vertex (see Fig. 1) to establish signals that can be 
related to semileptonic , B^, B^, and decays. The K^K~ vertex reconstruction 
does not use a constraint to the known mass [42], although |m(0) — 1.019| [GeV/c^] is 
required in order to select a pure sample of 4> candidates. 

3. Backgrounds to Lepton- Charm Signals 

Several backgrounds affect the semileptonic B signals. Some of these can be reduced 
by judicious signal selection, while some must be included in the modeling of the signal or 
treated as sources of systematic uncertainties. The simplest of these backgrounds to un- 
derstand are those events arising from combinatoric sources, which are generally estimated 
from the sidebands of the charm signal. In these backgrounds, random tracks are combined 
to form a charm signal which passes all charm selection requirements. This combinatoric 
background can most easily be reduced by selection requirements and modeled by the side- 
band events, which are expected to exhibit the same shape underneath the signal. A related, 
but more subtle type of background is that arising from the mis-identification of tracks in 
one charm decay arising from incorrect assignment of particle identifications in a real charm 
decay, resulting in "reflection" backgrounds. These backgrounds are often flat beneath the 
signal of interest, but occasionally they exhibit particular shapes that can affect the signal 
distribution non-uniformly. Some rcffcction backgrounds can be effectively reduced with 
particle idcntiflcation selections, such as the speciflc ionization of particles, dE/dx (see Sec- 
tion lie 4.) Other reflection backgrounds, which have non-uniform distribution in mass 
beneath the charm signal are included in the fit to the signal (see Section II C 5.) MC 
simulated data is used to determine the shape of these reflection backgrounds. 

The third type of background to the semileptonic signals arises from physical processes 
that produce a real lepton and charm hadron, but not through a decay directly to l~D. This 
includes processes which originate from the same B, such as S — > DD, where D — > l~X, 
and B t^Vt-D, where t' — > I'X. These "physics backgrounds" are included in the 
fit to the sample composition (see Section III). Other backgrounds include processes in 
which the lepton and charm hadron originate from separate hb and cc quark pairs, i.e. 
b DX, b £~^X, or c ^ DX, c — > £~X. The bb background gives a wrong sign (WS) 
lepton-charm combination, in which the charm and lepton have the same charge, while the cc 
background gives right sign (RS) lepton-charm combinations, in which the charm and lepton 
have opposite charge. All of these processes arc also possible with a real charm hadron and 
a false lepton. In the case of false leptons, both right sign and wrong sign lepton-charm 
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are expected to be present. Backgrounds which do not originate from the same B hadron 
are treated as a source of systematic uncertainty and described by the wrong sign lepton- 
charm events, which primarily describe false leptons (see Section VI A.) The cc background is 
assumed to be small for a charm decaying to a lepton with Pt > 4: GcV/c [43] and is ignored, 
while the bb background is implicitly included in the false lepton systematic uncertainty. 



4- Signal Optimization 

Requirements to further enhance the lepton-charm signal include px cuts on the p, K, 
and TT charm daughter tracks, and cuts on the invariant mass of the lepton-charm system, 
m{i~D) = \/{p^{i^) + p'^{D)Y, to limit feed-down from excited charm and lepton-charm 
combinations which do not originate from direct semileptonic bottom hadron decays. Re- 
quirements are also made on the probability of the charm and lepton-charm vertex fits. 

Since bottom hadrons are longer-lived, a powerful discriminant against these backgrounds 
is a cut on the proper time of l~D candidate. The decay distance of the B hadron is 
determined by defining a quantity, Lxy(PV i~D), which is the transverse decay distance 
of the lepton-charm combination from the primary interaction vertex (PV), projected on 
the i'D momentum direction. The missing neutrino produced in the semileptonic decay 
prevents precise knowledge of Pt{B) and thus of the proper decay time of the B candidate. 
Instead, a pseudo proper decay time is constructed as: 

ct*{rD) ^ L.y(PV ^ i-D) X (1) 

A ct*{£~D) > 200 yum cut is applied to guarantee a signal from long-hved bottom hadrons 
and to reduce signal contamination from false leptons and other processes that can contribute 
a lepton and a charm hadron from uncorrelated sources (see also Section II C 3.) This 
requirement also drastically reduces the combinatoric background of charm candidates with 
real leptons. A cut on the significance of the transverse decay distance of the charm meson, 
Lxy(PV — > -D)/cri^y(pv_>D), also reduces the light flavored hadron contamination in the 
signal. A cut on pt{D°) > 5 GeV/c is apphed to improve agreement between the £~D^ data 
and Monte Carlo simulation used in determining the efficiencies (see Sec. IV). The selected 
D*~^ candidates are a subset of the candidates. Instead of performing a vertex fit on the 
soft pion, 7r+, from the D*+ D^n+ decay, A tight Am{D*+,D^) = m(L'°7r+) - m(D°) G 
[0.1440,0.1475] GeV/c^ cut is used to select a very clean £~D*^ sample. This reduces the 
systematic uncertainty in the selection of the £~D*^ combination relative to a pair, 
since no additional vertex fit is performed. Consequently, the efficiency to detect the soft 
pion is better described by the simulation. Since the data agrees well with the simulation 
for tracks with pt greater than 400 MeV/c, as can be seen in Fig. 2, the soft pion efficiency 
is determined from the simulation. A tight Am{D*+,D'^) e [0.1440,0.1475] GeV/c^ cut is 
used to select a very clean £~D*^ sample. 
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FIG. 2: Comparisons between data and simulation of ptC^t^), the soft pion from the D*~^ decay, 
for the (a) fj.~D*~^ and (b) e^D*+ mode. 



In order to determine the final analysis selection, kinematic selection criteria are opti- 
mized with respect to the combinatoric background for each lepton-charm channel, with ad- 
ditional cuts designed to limit non-combinatoric background, such as the ct*{i~D) snadpriD) 
cuts, applied during the optimization. The figure of merit (FOM) used for optimization is 
S/ yS + B. The signal, 5, is taken from inclusive B — > i'u^DX and A° —>■ I'v^A.'^X Monte 
Carlo (see Sec. IV A). The background, -B, is taken from the sidebands of the charm signal. 
In order for the FOM to accurately reflect the significance of the signals in data, S is scaled 
to the expected data signal with a set of nominal cuts obtained by first optimizing each cut 
individually without applying any other cut. The cuts are then optimized a second time 
applying all optimal cuts from the prior optimization except the cut being optimized. After 
two or three successive iterations, a stable optimal cut point is reached for all cuts. 

A particle identification cut using dE/dx is found useful for reducing the combina- 
toric background in the A+ signal. The combinatoric background can be significantly 
reduced by correctly identifying the proton from the A^ pK~7[^ decay utilizing the 
specific energy loss of the proton track measured in the COT. A dE/dx likelihood ratio, 
CTZ, requirement is applied to the proton. The likelihood ratio is defined by the relation 
Cn{p) = C{p)/[C{p) + C{K) + C{7r) + C{e) + £(/!)], where C{i) oc exp{-Zf/{2a%)} and 
Zi = ln.[dE /dx^'^^^' /dE /dxY'^'^']- Figure 3 shows the resulting CTZ distributions for protons 
from the A'' pn^ decay and kaons and pions from the D*^ — > (— K~7r~^) vr"*" decay 
with the proton hypothesis applied. Muons are indistinguishable from pions, while electrons 
are well-separated from all of the other distributions, since their mass is so much lower than 
the mass of the other particles. A cut on CTZ{p) > 0.3, as determined from the control 
samples, is applied to reduce background while keeping the proton efficiency high. This cut 
primarily removes pions, since the dE/dx separation between protons and kaons is not as 
good. 

To cancel as many differences in signal reconstruction as possible, the selection criteria 
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FIG. 3: dE/dx CIZ distribution for protons from — > pvr and kaons and pions from D*+ 

{-^ K^TT^) 7r+ with the proton hypothesis apphed. Tracks with CTZ to the right of the dashed 
vertical hne are identified as protons. 

are kept as similar as is feasible across charm channels. The optimized cuts designed to limit 
both the combinatoric and some non-combinatoric backgrounds are unified to minimize the 
differences in selections between channels. However, some cuts, in which different optimal 
values are expected due to differences in the decay kinematics, are not forced to be similar. 
For example, the proper decay time of the meson and baryon differ by a factor of 
about five. The selection criteria applied to the lepton-charm decay signatures are listed 
in Table I. Additional selection requirements to reduce non-combinatoric backgrounds are 
discussed next. 



5. Reflection Backgrounds 

The selection criteria discussed above (see Sec. II C 4) optimize the signal sensitivity with 
respect to the combinatoric background. However, there are other non-combinatoric back- 
grounds that must be considered. This is partially achieved with the ct*{i~D) and Pt{D) 
cuts discussed previously. Another significant background arises from reflections, which oc- 
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TABLE I: Signal selection requirements. 



Selection cuts 












ct{D) [cm] G 


[-0.01,0.10] 


[-0.01,0.10] 


[-0.01,0.20] 


[-0.01,0.10] 


[-0.01,0.05] 


ct*{i-D) [cm] > 


0.02 


0.02 


0.02 


0.02 


0.02 


act*{e-D) [cm] < 


0.04 


0.04 


0.04 


0.04 


0.04 


m{e~D) [GeV/c2] G 


[2.4,5.1] 


[2.4,5.1] 


[2.4,5.1] 


[2.4,5.1] 


[3.4,5.5] 


Pt{D) [GeV/c] > 


5.0 


5.0 


N/A 


N/A 


N/A 


Pt{p) [GeV/c] > 


N/A 


N/A 


N/A 


N/A 


2.0 


Pt{K) [GeV/c] > 


0.6 


0.6 


0.6 


0.6 


0.6 




10 


10 


10 


10 


5 


vertex pToh.{£~ D) > 


10-7 


10-7 


10-7 


10-7 


10-4 


Lxy/cTLxyiD) > 


4.5 


4.5 


11 


5 


4.5 


Am{D*+,D^) [GeV/c2] g 


N/A 


[0.1440,0.1475] 


N/A 


N/A 


N/A 


Pt{tt*) [GeV/c] > 


N/A 


0.4 


N/A 


N/A 


N/A 


|m((/)) - 1.0191 [GeV/c2] < 


N/A 


N/A 


N/A 


0.0095 


N/A 


dE/dx Cn{p) > 


N/A 


N/A 


N/A 


N/A 


0.3 



cur when the particle identifications in charm decay are mis-assigned. For example, if the 
K'^ from a Df — > K'^K~7r'^ decay is assigned the pion mass, the n'^K'n'^ combination 
can contribute to the signal. Figure 4 shows the shapes determined from MC for refiec- 
tions from (a) (b) (c) D'^, (d) Df, and (e) A+ decays when these decay channels 
are reconstructed as a different charm mode. The shapes are normalized to their expected 
contributions, e.g. assuming fu '■ fa '■ fs '■ /a;, — 0.4 : 0.4 : 0.1 : 0.1, where these numerical 
values are for illustrative purposes only. The D'^ — > K^K~t:'^ decay is the most significant 
reflection background below the signal, shown in Fig. 4(c). This reflection is particularly 
problematic because the reflection begins just underneath the real signal. Potential 
p-TT mis-identification is a significant consideration in the signal, shown in Fig. 4(e). The 
K~Ti^Ti^ decay significantly contributes to the background beneath the A+ signal, 
although its contribution is fiat underneath the signal. 

The shape of the — > K^K~'k^ refiection background beneath the K~n~^7i~^ sig- 

nal is determined from a Monte Carlo simulation (see Section IV A) study, in which semilep- 
tonic Bg i~PiD'^X decays are generated. In these MC events — > K~7r~^7r~^ candi- 
dates are then reconstructed. The resulting K'tt'^tt'^ invariant mass distribution is shown 
in Fig. 5. The normalization of the — > K'^K^n'^ reflection shape in the flt to the 
signal is determined by reconstructing a Df — > 0(— > K~K~^) n'^ signal from the wide 
signal window, 1.78 GeV/c^ < m{D+ K-7r+7r+) < 1.95 GeV/c^ shown in Fig. 6. A 
mass cut of \m{K'^K~) — 1.019 GeV/c^| < 0.0095 GeV/c^, designed to reduce background 
to the Df signal, is applied to the (f) — > K'^K~ decay. Monte Carlo simulation is then 
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FIG. 4: Monte Carlo simulation reflection shapes for (a) (b) (c) (d) D+, and (e) A+. 

The shapes are normalized to their expected contributions, assuming fu'-fd-fs- = 0.4 : 0.4 : 
0.1 : 0.1, used for illustrative purposes only. 



used to measure the efficiency of the —>■ (p-n^ decay relative to the inclusive set of 
K^K~Ti^ decays that contribute to the reflection. The converse K~ti^tt^ re- 

flection in the K'^K^tt^ signal is negligible due to the mass cut applied to the 

K^K~ invariant mass. 
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In a manner completely analogous to the way the signal yield, A^data(-D^ 07r+), 
is determined in data, the candidates decaying to the 07r"*" and K^K^n^ states, 
Nmc{D^ 07?+) and Nmc{.D^ K^K^n^), respectively, are determined from the Monte 
Carlo simulation. The number of mesons expected to contribute to the signal can 
then be calculated by evaluating 



iVdata(/^+ ^ K^K- 



TX 



iVdata(/^+ ^ 07r+) 



(2) 
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where 



The numbers of candidates that contribute to the lepton-charm samples in the 
wide mass window around the signal are Ne{Df — > K^K^ti^) = 1580 ± 130 and 
N^{Df K+K-TT+) = 2570±210. The normahzation of the Df reflection in the D+ signal 
is later constrained to the predicted number of Df reflection events in the flt to the signal 
(see Sec. II D). 

Since the D'^ and Df reflections in the A+ signal are relatively flat under the signal 
region, sideband subtraction is expected to remove the effect of the and Df reflections 
on the A+ signal distributions within statistical uncertainty. Correspondingly, the event 
count obtained by fitting the A+ signal is not expected to be signiflcantly influenced by 
the presence of these backgrounds. Additionally, the dE/dx cut applied to the proton 
(discussed in the previous section) reduces contamination from pions, which contribute to 
the —>■ i^"7r+7r+ and D*+ [i^"7r+]7r+ reflections. 



D. Signal Yields 

The m{K^ 71^71^), m{K~pi'^), m{K'^K^pi~^), and m{pK^pi^) mass spectra are fit to 
determine the number of lepton-charm events for the i^D^, i^D^, i^D*^, i^D^, and ^^A+ 
samples. The invariant mass distributions of the charm signals are shown in Fig. 7 for 
the /x+SVT data and in Fig. 8 for the e+SVT data with all lepton, charm, and lepton- 
charm selection criteria apphed. The distributions are flt with a double Gaussian and linear 
background shape. The reflection of Df decays into the flnal state is included in the 
flt to the signal. The normahzation of the Df reflection is constrained to the predicted 
number of Df reflection events as described above. In order to keep the broad Gaussian 
and reflection shapes reasonably independent, the double Gaussian means and widths for the 
D'^ are determined before the reflection shape is added to the flt. When the combined flt is 
performed, the parameters of the double Gaussian are constrained within their uncertainties. 
The flts to the D^, D*^, D^, D'^, and A+ charm signals for right sign lepton-charm pairs 
are shown in Fig. 7 for the /x-l-SVT data and in Fig. 8 for the e-|-SVT data. The invariant 
mass distributions for wrong sign combinations of lepton-charm pairs, e.g. D^£^, are also 
included in Figs. 7 and 8, indicating no significant contributions of possible backgrounds, 
such as false leptons, to be present in the right sign signals (see also Sec. VIA). The fitted 
lepton-charm yields are listed in Table II. The reflection is not included in the yield, 
since the fit shape to the includes a separate shape for the D'^ refiection, as discussed 
in Section II C 5. The dE/dx cut flattens the background and reduces its overall level by a 
factor of flve, while it reduces the signal by ~ 35% in the /i-|-SVT data and ~ 28% in the 
e-FSVT data as can be seen in Fig. 7(e)-(f) and Fig. 8(e)-(f). 
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FIG. 7: ;U+SVT right sign (RS) (points with error bars) and wrong sign (WS) (histogram) invariant 
mass distribution of (a) D^, (b) Z)*+, (c) , (d) Df, (e) A+ with ah cuts apphed and (f) without 
the dE/dx cut apphed. The fit parameterizations described in the text are overlaid. 

III. SAMPLE COMPOSITION DETERMINATION PROCEDURE 

This measurement uses flavor SU(3) symmetry to describe the branching fractions of 
semileptonic B meson decays; therefore, the partial widths of the semileptonic decays of 
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FIG. 8: e+SVT right sign (RS) (points with error bars) and wrong sign (WS) (histogram) invariant 
mass distribution of (a) D^, (b) D*~^, (c) L*^, (d) Df, (e) with all cuts applied and (f) without 
the dE/dx cut applied. The fit parameterizations described in the text are overlaid. 

B mesons are chosen to be equal, namely 

r(5° ^ rueX) = T{B- e-UiX) = r(5° ^ r^eX), (4) 
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TABLE II: Fitted signal yields for lepton-charm final states in 360 pb 







e±J3 V 1 






/i±o V 1 




Decay 


Yield 


FOM 


Fit Prob. [%] 


Yield 


FOM 


Fit Prob. [%] 


r^o 


16, 939 ± 160 


122 


64.4 


29, 909 ± 224 


159 


12.5 


£-D*+ 


2, 998 ±56 


54.1 


1.27 


5, 492 ± 77 


73.3 


1.14 


£-D+ 


10, 779 ± 149 


90.2 


9.43 


20, 236 ±216 


114 


50.7 


l-Dt 


1,012 ±44 


27.3 


7.84 


2, 069 ± 84 


36.6 


30.2 




1, 755 ± 106 


32.8 


33.9 


2, 984 ± 130 


40.9 


40.9 



where 

V{B ^ ru,X) = :^^13{B ^ i-i^eX). (5) 

This assumption is referred to as the spectator model, which also implies that the partial 
widths of the semileptonic bottom hadron decays into the pseudoscalar, vector, or higher 
excited D states are expected to be equal, 

r(BO ^ i-iyeD+) = T{B- ^ rPeD^) = r(50 ^ rhDj) = T{B ^ ri/^D), (6) 

with similar relations holding for D* and D** decays. The additional constraint that 

r{B rueD) + r*{B^ r^eD*) + r**{B rviD**) = v{b rviX), (7) 

is also applied to the partial widths. This constraint includes non-resonant decays and 
h ^ u transitions in addition to actual D** decays in the D** partial width. Since excited 
A° semileptonic decays are not necessarily well-described by the spectator model, fixed 
branching fractions are used to describe those decays [44] (see Table VI) . 

A simplified example illustrating the procedure used to extract the sample composition 
follows. Assuming that the only source of £~D'^ combinations is from the direct decay of a 
neutral B meson, such as B^ — > l~ViD'^, the number of reconstructed I'D^ events can be 
expressed as 

N{e-D+) = N{B°) X B{B° rufD+) x B{D+ K-ti+ti+) 

= [N{B) ■ fd] X [r(5°) ■ r(5° ^ rz/^D+)] X B{D+ i^-7r+7r+) 

X £(S° ^ ri/^D+, D+ K-7r+7r+). (8) 

The number of reconstructed i~D'^ combinations, N{i~D^), can be related to the number 
of B^ mesons, iV(i?°), produced in the fragmentation process by the branching fraction of 
the S° — > i'PiD'^ decay, the branching fraction of the K'n'^n'^ charm decay, and 
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FIG. 9: Illustration of cross talk between S-meson species. 



the detection and reconstruction efficiencies for the entire decay chain. N{B'^) and N{B) 
represent the number of and generic bottom hadrons produced, respectively. 

However, as sketched in Fig. 9, cross-talk between the various bottom hadron species via 
the excited charm states necessitates a sample composition parameterization to relate the 
lepton-charm signals to the parent bottom hadrons. Introducing the relative fragmentation 
fractions, the sample composition for the B^ can be written (for illustrative purposes) as: 



X 



.BO-^D+X 



-ueD+X)e{B- 



yiD+X) 



B-^D+X 
fs 



fu + fd 



^ •''^^ BO^D+X 



■ (9) 



The sample composition parameterization requires knowledge of the branching fractions 
of the charm hadrons, which are determined from the world average values as compiled by 
the Particle Data Group (PDG) [42], when available, and from theoretical predictions and 
symmetry principles [45], when not available in the PDG. The parameterization also re- 
quires knowledge of the efficiency of reconstructing the decay, which is primarily determined 
from Monte Carlo simulations. In addition to the primary decays that contribute to the 
semileptonic signal, indirect semileptonic decays {e.g. B DDX, D i~X) contributing 
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to the l~D final state are also included in the parameterization of the sample composition. 
All of the decays considered in the bottom hadron sample composition procedure are listed 
in Table III (see also Fig. 9). 

The bottom meson branching fractions are included via the partial widths, listed in Ta- 
ble IV, and adjusted by the lifetime of the respective bottom hadron, given in Table V. The 
branching fractions used for the A° semileptonic decays are estimated from measurements of 
the branching fractions made in other CDF measurements [44], as shown in Table VI. The 
ground state charm branching fractions used in this measurement are listed in Table VII. 

IV. EFFICIENCIES 

Since the bottom hadron fragmentation fractions are measured relative to each other, 
most efficiencies in the measurements are expected to cancel. Many of the remaining relative 
efficiencies are determined from Monte Carlo simulated data in which the trigger and all 
detector calibrations are configured just as they are determined for a given run in real 
data. Comparisons between the data and the inclusive simulation samples, discussed in 
the subsequent section, validate the use of simulation to estimate the relative efficiencies 
between B semileptonic decays and the lepton-charm signals. A few absolute efficiencies, 
such as the different XFT trigger efficiencies for K, tt, and p are not properly described 
in the Monte Carlo simulation. These efficiencies must be determined from data and are 
discussed in Section IV B. 

A. Monte Carlo Simulation 

Monte Carlo simulation is used at various points throughout the measurement. Although 
the simulation utilized in the sample composition process is generated both for exclusive bot- 
tom hadron decays and inclusive B semileptonic decays, all of the Monte Carlo simulation 
samples used in the measurement have the same parameters for generation. The simu- 
lated events are passed through the GEANT-based [46, 47] CDF II detector simulation [48]. 
GEANT 3 simulates the passage of the long-lived particles through the material of the de- 
tector and includes multiple scattering effects. All simulated samples are generated with a 
"realistic" , rather than parametric, simulation. A tuned magnetic field and Geant material 
description [49] are applied in order to correct for regions of the detector where the material 
is under-represented in the simulation. 

A single bottom hadron is generated according to an input transverse momentum and 
rapidity spectrum, which have been determined from data. The pt spectrum, obtained in 
the inclusive J/%1) cross-section measurement [20], is used as the input B meson spectrum. 
However, it appears that the momentum distributions of the A° decay products, in particular 
the ^~A+ momentum spectrum, are not well described using the same spectrum as is used for 
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TABLE III: Bottom hadron semileptonic sample composition used in the measurement. NR refers 
to non-resonant decays. 

go B Ag 

£-iyD+ i-uD" ri>L>+ rpA+ 

e^i?D*+ e-uD*^ i-i?D*+ £-i>Ac(2593)+ 

^D°Tr+ ^Dfl ^ Sc(2455)++7r- 

r>+7r77 A+7r+ 

^ Sc(2455)07r+ 

rz^D? ^-i>L>+i(2460) -^A+TT- 

^L>*07r+ ^ZJ^V ^-C)^(]"7r° ^ Se(2455)+7rO 

^L'*+7r° -^D*+Tr- ^ D+j ^ A+Tr+vr" 

D+TT^f-f D+TT^/j A+7 



£-pD^o i-uD*+ (2317) i-uAc{2625)+ 

-^D%+ ^L>V ^-D+ttO -.A+tt+tt- 



A+7 



^D^TT^j -^L'V/7 ^D^Tr+ ->A+7r+ 

^ L>*+7rO ^ D*+Tr- D+7r°/7 

L>07r+ D'^7r+ D*^K+ £-9^^2455)^7?+ 

L»+7r0/7 D+TT^/j -^DV/7 ^A+tt- 

£-i?D'+{2573) ^-pSc(2455)+7r° 

^L>*07r+ ^D^OttO ^D*+K^ ^A+TfO 

^L>*+7rO ^L'*+7r- L'+7r77 £-vA+fQ 

L>+7r77 D+tt'^/j -^D07r°/7 i>A+7r+7r-(NR) 

^£>+7r- ^D+TT- ^D^K+ ri^A+7rV(NR) 

£-i?D*+'iT^ (NR) £-i?D*+Tr-{NR) ^-pD*+7rO(NR) 

-> D07r+ ^ D07r+ ^ i:»+7 

Li+7r77 D+TT^/j 

^-i?D+7rO(NR) 

rpL'*°7r+(NR) rpD*°7r°(NR) 

^ L>°7r77 ^ L»07r77 

i'-i^D+vr^NR) ^-z>i:>+7r-(NR) 

£-z?i:)°7r+(NR) ^-i?D07r°(NR) 

£){*)£){*) X D^*^D^*^K D^*^D^*^K 

£)(*)+ 

^(*)^(*)^ D^dWa: r-jyA+ 

r-i/Ac(2593)+ 

r-z/D+W'(**) r-i/i:»°W'(**) T-i/i:>^^*^'^**^ t-i^Ac(2625)+ 
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TABLE IV: Partial widths of the B mesons used (from Ref. [42]). 



B Decay 


Partial Width [ps"^] 


T{B e-i^iD) 


0.0134 ± 0.0009 


T*{B e-i?iD*) 


0.0372 ± 0.0017 


T**{B ^ £-i?eD**) 


0.0141 ± 0.0010 


TABLE V: Bottom hadron hfetimes used 


in the measurement (from Ref. [42]). 


B Lifetimes 


[ps] 


rm 


1.536 ±0.014 


t{B-) 


1.671 ±0.018 


r{B-)/rm 


1.086 ±0.017 


rm 


1.461 ± 0.057 


r{K) 


1.229 ± 0.080 



the B mesons, as can be seen in Figure 10(a). The Monte Carlo simulation generated /i"A+ 
transverse momentum spectrum is observed to be harder than the measured semileptonic 
spectrum. This indicates a potential difference in the momentum dependence of b baryon and 
B meson fragmentation processes. Consequently, instead of using the spectrum used for the 
mesons, apT spectrum derived from the semileptonic £~A+ data is used in the measurement. 
This tuned spectrum [43] shows good agreement between data and Monte Carlo simulation, 
as shown in Fig. 10(b). The tuned ^~A+ spectrum is obtained by re- weighting the bottom 



TABLE VI: A° branching fractions B used in the measurement (from Ref. [44]). The lack of a 
quoted uncertainty indicates an assumption made for B. NR refers to non-resonant decays. 



A|] Decay 


B 




(9.2±2.1)% 


ri>^A+ 


(6.54±0.22)% 


£-i>^Ac(2593)+ 


(3.07±1.02)xl0 


£-i>^Ac(2625)+ 


(5.14±0.99)xl0 


£-i>fSc(2455)++7r- 


(2.7±1.0)xl0-^ 


rz^^Sc(2455)+7r° 


(2.7±1.0)xl0-^ 


ri^^Sc(2455)07r+ 


(2.7±1.0)xl0-3 




2.6x10-3 


i-VlkcTT+TT- (NR) 


5.2x10-3 


^-z^£Ac7r°7r° (NR) 


2.6x10-3 
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TABLE VII: Ground state charm branching fractions B used in the measurement (from Ref. [42]). 



Charm Decay B [%] 

D+^K-TT+TT+ 8.8 ±0.6 

L»o K-TT+ 3.80 ± 0.09 

D*+ ^ D°TT+ 67.7 ±0.5 

D+^(I)TT+ 3.6 ±0.9 

(p^K+K- 49.1 ±0.6 

A+^pK-TT+ 5.0 ±1.3 
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FIG. 10: Comparisons of /i~A^ transverse momentum spectrum between data and Monte Carlo 
simulation generated according to the pT spectrum inferred from (a) the inclusive J/^l' cross-section 
measurement and (b) the tuned semileptonic spectrum. The corresponding bottom plots show 
the ratio of data over MC with fits of a constant (dotted) and straight line (solid) overlaid. 



hadron pt spectrum measured from the inclusive J/tp cross-section measurement [20]. The 
re-weighting function is determined from the disagreement between the ^~A+ data and 
the generated spectrum, which is fit to a first order polynomial, w = b + m ■ pt, where 
b = 1.43 ± 0.08 and m = —0.026 ± 0.007 are the values of the fit averaged between the 
e±SVT and ;U±SVT data. The tuned spectrum is then varied by ±2 a of the uncertainty on 
the slope (am) to bound the uncertainty on this spectrum (see Sec. VII). All of the Monte 
Carlo simulation events are generated with an input pt threshold of pt(-S) > 5 GeV/c and 

\vm < 1.1. 

After the bottom hadron is generated, it is then decayed using the EvtGen decay pack- 
age [45], which decays the particles according to a user-specified decay chain and theoretical 
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decay models. The B meson form factors used in the Monte Carlo simulation for this anal- 
ysis are taken from various models based on heavy quark effective theory (HQET) [50, 51]. 
The ISGW2 [50] model implemented in EvtGen governs the B meson semileptonic de- 
cays to the ground state and doubly excited charm mesons, while the HQET decay model, 
implemented in EvtGen, is used for the B meson semileptonic decays to excited charm 
states. Non-resonant D** meson decays are described by the model developed by Goity and 
Roberts [51]. The Aj) baryon semileptonic decay model is newly implemented [43] into the 
EvtGen package for this measurement. The baryon form factors for the primary semilep- 
tonic £~i/£A+(*'**^ decays are taken from constituent quark model calculations made 
by Pervin et al. [52]. These results agree with the large Nc predictions by Leibovich and 
Stewart [53] to order ©(I/ttiq). Non-resonant A° decays, which are expected to contribute 
comparatively little to the total A° semileptonic width, are described by a phase space decay 
model. 

1. Data - Simulation Comparison 

Four inclusive Monte Carlo simulation samples, S° — > ^~P^D°'+X, B~ I'V^D^^'^X, 
Bg — > £~9iDfX, and A° — > £~P^A+X are generated to vahdate the use of simulation to 
determine the kinematic efficiencies of the bottom hadron semileptonic decays used in the 
measurement. The agreement between the data and the Monte Carlo should not be very 
sensitive to variations in the D** branching fractions between the default EvtGen table 
and the one to be later determined in the fit for the fragmentation fractions. The agreement 
between data and Monte Carlo simulation is checked for quantities used in the signal selection 
(listed in Table 1). In general, the agreement between data and simulation is good in both the 
/i+SVT and e+SVT data. A typical example of comparisons between data and Monte Carlo 
in the ;U+SVT sample is shown in Fig. 11 for (a) ct* {n^ D^) , (b) act*{lJ'^D^), (c) prilJ^^D^), 
and (d) m{n~D~^). A complete set of comparisons between data and MC can be found in 
Ref. [43]. The area of the simulation distribution is normalized to the corresponding area 
of the data distribution for this comparison. The quality of the comparisons are quantified 
by fitting the ratio of the data to the simulation by both a first order polynomial and a 
constant. The former indicates potential biases between the two distributions {i.e. whether 
the simulation distribution is too hard or soft relative to the data), while the latter gives 
a measure of overall agreement between the distributions. No significant disagreement, 
determined from the fit to a constant line, is observed between data and simulation in the 
quantities used for the signal selection. 
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FIG. 11: Data-MC comparisons for fi-D+ of (a) ct*{^-D+), (b) act*{f^-D+), (c) pt{^i'D+), and 
(d) m{ix-D+). 

B. Relative Efficiency Determination from Data 



Many efficiencies in the measurement of the relative fragmentation fractions are expected 
to cancel, and many of the remaining relative efficiencies are determined from the Monte 
Carlo simulation. The few efficiencies that are not well described by the simulation are 
determined from the data and discussed next. 



1. XFT Efficiency 

Differences in the XFT efficiencies of kaons, pions, and protons are expected due to the 
stringent hit requirement placed on COT tracks by the XFT trigger. Since the dE/dx of 
kaons and protons is lower than the di?/dx of pions for pt > 2 GeV/c, the COT hit require- 
ment leads to a lower efficiency for kaons and protons relative to the pion XFT efficiency. 
These efficiencies are difficult to describe in the Monte Carlo simulation due to varying COT 
operating conditions during the data-taking period of this measurement. Therefore, they 
are derived from the data. The SVT efficiencies, which contribute to the triggers used in 
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TABLE VIII: Parameterizations for the XFT Monte Carlo simulation corrections. 







K 




TT 


Data period 


ao 


Ol 


ao 




Fcb"02-Oct'02 


0.9931±0.05 


-0.0725±0.02 


0.9772±0.03 


0.00968±0.01 


Oct'02-Jan'03 


0.9584±0.02 


-0.1952±0.007 


1.0016±0.01 


-0.1501±0.005 


Jan'03-Jun'03 


0.9359±0.02 


-0.1919±0.007 


0.9851±0.01 


-0.1341±0.004 


Jun'OS-Sep'OS 


0.9282±0.01 


-0.1897±0.005 


0.9921±0.008 


-0.1776±0.004 


Sep'03-Aug'04 


0.9643±0.01 


-0.0907±0.004 


0.9931±0.007 


-0.0678±0.003 








proton 






ao 


ai 




as 


Feb'02-Jun'03 


1.063±0.090 


-1.326±0.963 


3.198±3.218 


-2.203±3.391 



this measurement, depend directly on the XFT efficiencies. The species dependence of the 
SVT efficiencies originate entirely from the XFT, since the energy loss between K, ir, and 
p are negligible in the silicon detector relative to the drift chamber. The differences in the 
SVT efficiencies between reconstructed lepton-charm channels are therefore described by the 
dependence of XFT efficiencies on particle species. 

The XFT efficiencies for K and tt are measured by reconstructing the D'^ — > K'n'^n'^ 
decay mode in the two-track trigger (TTT) data sample, where two tracks are required to 
match to the SVT trigger and no lepton requirement is made. Two of the final state K 
or TT tracks are matched to the SVT tracks. The track that is not matched to an SVT 
track is treated as the unbiased track, which is then examined to determine whether it could 
have fired the XFT trigger. Tracks that could have passed the XFT trigger are included in 
the numerator of the efficiency, while all unbiased tracks are included in the denominator. 
A similar procedure is carried out for the proton XFT efficiency, using A" — >• pTT events 
reconstructed in data collected with the TTT, where two other tracks in the event are 
required to have ffied the SVT trigger. These efficiencies, binned in time to span the data 
set used, are shown in Fig. 12 for K and tt. The ratio of the K and tt efficiencies determined in 
the data relative to those determined in the corresponding Monte Carlo simulation are shown 
in Fig. 13. These corrections are parametrized by linear functions of the form ao + cli/pt-, 
and the obtained fit parameters are listed in Table VIII. Details of the fits for the proton 
XFT efficiency can be found in Ref. [43]. 

2. dE/dx Efficiency 

The dE/dx efficiency of the CTZ cut apphed to the proton, discussed previously in Sec- 
tion lie 4, is also evaluated from data and the Monte Carlo simulation is adjusted accord- 
ingly. The A° — > p7r~ control sample is used to evaluate the efficiency of the CTZ{p) > 0.3 cut 
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FIG. 12: XFT efficiency exFT as a function of pj.^ for pions (upper curve) and kaons (lower curve) 
for various run ranges indicated in distributions (a) througli (e). 

applied to the proton from the decay. The dE/dx efficiency is obtained by dividing the 
number of protons that pass the dE/dx CTZ{p) cut by all protons in bins of proton transverse 
momentum. The shape of the efficiency is parameterized by two functional forms: a third 
order polynomial plus a constant and using only a constant. Both parameterizations fit the 
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FIG. 13: Ratio exFT(DATA)/exFT(MC) as a function of pj, for pions (upper line) and kaons 
(lower line) for various run ranges indicated in distributions (a) through (e). 

data well; the former is used as the default parameterization in the measurement, while the 
latter is used in the evaluation of systematic uncertainties (see Sec. VI). 
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3. Single Track Efficiency 

The difference in efficiency between the two traclc charm topology in the — * K^tt~^ 
decay and a three track topology such as the — > K~7r~^7r~^ decay arises from the ef- 
ficiency of reconstructing an additional single track;. Since the Monte Carlo simulation 
contains only the decay products of the generated bottom hadron and no additional tracks 
from the fragmentation process or underlying event, the efficiency of reconstructing a track 
in a simulation event is different from that in data. Thus the single track efficiency is 
determined from the data relative to the simulation. In order to evaluate this efficiency, 
the lepton plus four track state with D° — > K'n^n'Tr'^ , is reconstructed and nor- 

mahzed to the ijl~D^ decay with — > K~7r'^ in both the data and the simulation. This 
represents the square of the efficiency, e^^^, to find a single track in the data relative to 
the simulation, assuming that the reconstruction of the third and fourth tracks in the 
D° — > ir~7r"'"7r~7r"'" decay are uncorrelated. This assumption will be treated as a source 
of systematic uncertainty, as discussed in Sec. VI. This procedure yields the single track 
efficiency in data relative to the same efficiency in MC, etrk = 87.8 ± 0.8 (stat) {B)%, 
where the systematic uncertainty is dominated by the knowledge of the relative branching 
fraction K' ti+ ti' ti+) / BTZ{D^ K~Ti+) = 2.10 ± 0.03 ± 0.06 [54], which is used 

to adjust the generated Monte Carlo simulation samples for both decay modes. 

4- Total Relative Efficiency 

The relative efficiencies Si-d that are included in the sample composition are the product 
of the acceptance £accep, lepton plus SVT trigger efficiency ^trigg, analysis efficiency ^an, and 
adjusted XFT efficiency Ep/K/n trigg, 

accep 

(MC) X £trigg(MC) X £an(MC) 

X Sp/K/n trigg(data/MC). (10) 

The relative efficiencies arc similar for the other lepton-charm signals. An extra factor of 
e^^^]^(data/MC) is needed for ei.ci(^~-D°) and e^cii^^ D*^) to adjust the two track charm 
topology relative to the three track charm states, while the i'^A^ relative efficiency requires 
an additional efficiency correction £dE/da;(data) for the dE/dx cut imposed on the proton 
candidate. 

Monte Carlo yields in each channel are determined by fitting the simulation signal 
with a double Gaussian, analogous to the fits to the data. A single Gaussian is used 
to obtain the yield in the indirect lepton-charm decays, where the lepton originates from 
B T'DX, T~ £~X, or another charm {e.g. B DDX, D £~X), because the 
yields are generally quite low in these channels and are poorly described by a double Gaus- 
sian. The total efficiencies derived for l~D^, l~D*'^ , and l~D'^ channels are listed 
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in Tables IX- XII for each B mode considered in the sample composition. The efficiencies 
are calculated relative to a reference channel where the yield out of 10^ generated events is 
listed for each corresponding reference channel. All of the decays shown in the tables are 
generated separately and for each decay lO''' events are generated. The relative efficiencies 
for the primary semileptonic decays do not include the branching fractions. The efficiencies 
for the "physics backgrounds" include the semileptonic branching fractions of the D or r, 
as well as the B DDX sample composition. Since a fixed sample composition is used 
for the — > i~j?£A'^X decays, the efficiencies quoted in Table XIII have the excited charm 
baryon branching fractions from Table VI applied. 



V. FIT OF RELATIVE FRAGMENTATION FRACTIONS 

A x^-fit is used to extract the fragmentation fractions from the sample composition of the 
semileptonic bottom hadron decays reconstructed. The measured yields in the five lepton- 
charm signals are fit to the yields predicted by the sample composition procedure, and the 
decay rates (r^*'**)) of the B meson to the ground and excited states are constrained within 
their errors. The x^, which is minimized, is 



^2 I -^(^ -^j ) measured N(^i -Dj ) predicted | 

\ (ri^Ni^i. -Dj ) measured ) / 



- J- PDG \ , / ^ ^ ^ PDG \ I / ^ ^ ^ PDG 



+ +[- — ^ +{- — ^ , (11) 



where A^predicted is determined from the sample composition process and Admeasured is obtained 
from the data (see Tab. II). In order to fit in terms of better measured quantities, the 
number of predicted lepton-charm events are expressed in terms of N[B^), which is an 
overall normalization in the fit and not indicative of the physical number of B^ mesons in 
the data, and in terms of the B meson lifetimes relative to the lifetime of the B^. The 
predicted number of lepton-charm events used in the fit for the B mesons is expressed as: 

Ad(rA)predicted = V Ar(S°) |^ X T (5°) 
j=a,u,s 

XT{Bj)/T{B^) 5^ Tfc X 5^ B{D,km ^ A) B{Di) Si^km 
k m 

fj f I la 

j=d,u,s fu + fd\ fd, 

xt{Bj)/t{B'')J2^, X J2^iDjkm ^ Di)B{Di)eijkm, (12) 

k m 

where A = D+, D^, D+, Tk = T, T*, T**, and the sum over m applies to F** if there is 
more than one D** state or non- resonant decay that can contribute to the final state D^. In 



= E Nmy^(l + ^)xriB' 
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TABLE IX: Efficiencies from I Monte Carlo simulation. For display purpose, the efficiencies 



are given relative to the ground state charm mode with the yield of the semileptonic decay into 


the ground state charm mode included in 


parentheses. 




Decays 


e+o V 1 


/i+O V i 


Jj — > I V D 


i.Uoy • iU (iU,oyU±iiUJ 


1 Qnv 1 n^3 /■ 1 Q nvnj-i on^ 
i.oU/ • iU (^io,U ( UitizUj 


ri n*+ ( r)+^o /,,\ 

t UU yU TT / 7j 


i.UDiitU.UiO 


1 1 nn-un ni q 


1)— r, n+ /'^O n*+^ (' n+„^0 /^,\ 
I VU^ [TT JJ ^ )[JJ^Tr / J) 


U.D/Z±U.Uii 


n ^c^n_i_n ni n 


o—r, n*+/' n+^0\ 
t I'J-yQ \J-> TT ) 


U.DZOltU.UiU 


n vi o_i_n ni n 


B— r, n'+ /'^O n*+A f /„a 
I vD^ [ttU ^ )[JJ^7r /^) 


n ^on_i_n ni n 


n T/IQ-Ln nin 


t VL)2 [TT L) )\L) ' TT /7j 


n ^7Q_Ln ni 1 


n vc;Q_i_n nin 


t VJ-J2 \ ^ ) 


u.oyu±u.uii 


n 7SQ-4-n ni 1 


t VU TT yU TT /^) yi\ix) 


n /iQc_i_n fine 


n fiQC-un nno 






n 7(\AA-c\ m n 

u. / u^mu. uiu 


n(*) r)(*) i>r^+ y"\ 




n nn/i/i-un nnnQ 


r)(*)+ n{*)~ ( f r)+ y\ 


n nnQ9-i-n nnn/L 


n ni fio-i-n nnni=; 


n(*) n(*) YfeTi+ y\ 


n nno7-i-n nnno 
u.uuz < ±u.uuuz 


n nnfiO-un nnnQ 


T—D r)(*)**) (fr)+ Y\ 


n noi 9-i-n nnnc; 


n n989-i-n nnn? 


D— V if— nO/"—— 7~)*+\/ 7~)+_-0 /^A 

r> — uJJiyTT U ^ )\U^TT /^) 


n ^ctc_Ln ni 1 


n Tct/i-Ln ni 1 


I VU^ yU TT ) 


n fioo-un ni n 


n 707_i_n ni n 
u. / z < ±U.U1U 


t VU-y \TT JJ jyJJ TT /7j 


n ^Ti -Ln ni 1 


n 7ccQ_Ln ni 1 




n K/ifi-un ni 1 


n 7c^o-i-n ni 1 




n fififi-i-n ni i 


n 7f;Q-i-n ni 1 
u. / uomu.ui -L 


f~ur)*+ -TT- ( D+ir'^ /'^^ f'NR'l 
•C u ly n y^ It) V / 


n 401 _Ln nin 


n «9«_|_n nnq 




n KQ/i-i-n noQ 


n 7fi74-n ni n 


n(*) r){*) T<"{('D~^ 


n nnnc;-(-n nnni 


n nni "^-t-n nnni 


nV*) T^(*) Y ( [) T\-\- Y\ 

Dg L)^ 'A ylL) A j 


n nnno_i_n nnni 
U.UUUo±U.U(JUi 


n nnoo_i_n nnno 
U.UUoo±U.U(JUz 




0.0031±0.0002 


0.0045±0.0002 


ri>L'l|(2535)(i^0L>*+)(i?+7rV7) 


0.612±0.010 


0.695±0.010 


rPi^:+ (2573) (i^0i?*+)(L'+7r77) 


0.575±0.009 


0.665±0.009 


rpz);^(2573)(i:)+-ft:°) 


0.592±0.010 


0.691±0.009 




0.0024±0.0002 


0.0073±0.0003 


D^*^D'<*'^X{ID+X) 


O.OOllzbO.OOOl 


0.0037±0.0002 



that case the sum is weighted according to the branching fractions BTZ{Djkjn — Di) of the 
various contributing D** states. The efficiencies Eijkm refer to the corresponding absolute 
efficiencies of events obtained in a particular final state normalized to the generated number 
of MC events with Pt{B) > 5 GeV/c. These efficiencies are detailed in Tables IX through 
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TABLE X: Efficiencies in £ D Monte Carlo simulation. 



Decays 




e+SVT 


//+SVT 


B° 






0.970±0.012 


1.028±0.012 




e-i?D+{Tr°D*+){D^Tr+) 




0.616±0.008 


0.702±0.009 




£-i?D+{7r+D*^){D°7T°/-f) 




0.605±0.009 


0.682±0.009 




£-pZ)o+(D°7r+) 




0.577±0.008 


0.679±0.009 




i-i?D'+{Tr^D*+){D^7r+) 




0.640±0.009 


0.700±0.009 




ri^r>i+(7r+L>*0)(L»07r77) 




0.557±0.008 


0.631±0.008 




i~i^D*+{Tr^D*+){D^7T+) 




0.613±0.009 


0.707±0.009 








0.622±0.009 


0.696±0.009 








0.640±0.009 


0.745±0.010 




t vD*^Tr(lT^T\^) (NR) 




0.461±0.007 


0.562±0.008 




^-i/D*°7r+(L>°7r°/7) (NR) 




0.451±0.007 


0.578±0.008 




l-vD^-K^ (NR) 




0.518±0.008 


0.698±0.009 








0.0024±0.0002 


0.0084±0.0003 




£){*)+ Di*)- (^£1)0 X) 




0.0033±0.0002 


0.0074±0.0003 








0.0026±0.0001 


0.007±0.0002 




T-uD^*'**\iD^X) 




0.0225±0.0005 


0.0311±0.0007 








1.376 ■ 10-3 (13,760±120) 


1.535 ■ 10-3 (15,350±130) 




£-uD*^{D^TT^/-f) 




1.024±0.012 


1.057±0.012 








0.643±0.009 


0.717±0.009 




l-uD^^{n^D*+){D^'K+) 




0.607±0.009 


0.710±0.009 








0.619±0.009 


0.716±0.009 




ri^£>;o(7rOD*0)(L>07rO/7) 




0.648±0.009 


0.742±0.009 




£-i^Df{TT-D*+){D^Tr+) 




0.609±0.009 


0.709±0.009 




£-i>L'|°(7r°L'*°)(L'°7r°/7) 




0.638±0.009 


0.726±0.009 




£ uDf{'K Z)*+)(D"7r+) 




0.598±0.009 


0.710±0.009 




£ 1/1^2 (1/ tt ) 




0.662±0.009 


0.757±0.009 




/? 7~i*4- — / nO +\ /'i\TT~) \ 

^ VU^TT {D^TT^) (iNRj 




0.456±0.007 


0.580±0.008 




t VV TT (i/ TT /7j (NR) 




0.474±0.007 


0.597±0.008 




£-i?D^TT^ (NR) 




0.565±0.008 


0.745±0.009 




£){*)i5{*)K(^i:)0j^) 




0.0033±0.0002 


0.0104±0.0004 








0.0044±0.0003 


0.0109±0.0003 




r"PL»(*'**)(^L>OX) 




0.0403±0.0007 


0.0523±0.0009 




^ £-pL»^|(2535)(i^0L>*+)(D07r+ 


) 


0.541±0.008 


0.638±0.008 




£- vDX (2535) (ii:+L»*0) (DOttO 


77) 


0.519±0.008 


0.633±0.008 




£-j>L';+(2573)(i^OL'*+)(L»o^- 




0.513±0.008 


0.614±0.008 




rPi:»;+ (2573) (i^+Z)*0)(L'07r77) 


0.501±0.008 


0.621±0.008 




£-i>£»;^(2573)(D°i^+) 




0.537±0.008 


0.639±0.008 




i:>W5Wi(:(£D0j^) 




0.0022±0.0002 


0.0072±0.0003 








0.0014±0.0001 


0.0056±0.0002 
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TABLE XI: Efficiencies in e'D*+ Monte Carlo simulation. 



Decays 


e+SVT 


Af+SVT 




0.888 • 10-3 (8,880±100) 


1.068 • 10-3 (10,680±100) 


i~uDf{7r^D*+){D\+) 


0.642±0.011 


0.696±0.010 


1-UD'+{-k^D*+){D^ti+) 


0.680±0.011 


0.691±0.011 


£-iyD;+{TT^D*+){D^7r+) 


0.642±0.011 


0.698±0.011 


e-i?D*+7T°{D°Tr+) (NR) 


0.484±0.009 


0.552±0.009 


D(*W(*^K{eD*+X){D°Tr+) 


0.0010±0.0003 


0.0028±0.0002 


D^*)+D^*)-{lD*+X){D^iT+) 


0.0033±0.0003 


0.0064±0.0004 


Di*^D(*^X{W*+X) (£>07r+ ) 


0.00105±0.00004 


0.0023±0.0002 


T-uD(*'**\eD*+X){D°TT+) 


0.018±0.003 


0.0235±0.0007 


B- i-uD^{TT-D*+){D^TT+) 


0.646±0.011 


0.696±0.011 


£-uD[^{Tr-D*+){D'^Tr+) 


0.640±0.011 


0.701±0.011 




0.620±0.011 


0.688±0.011 


i-vD*+iT-{D^TT+) (NR) 


0.466±0.009 


0.577±0.009 


D^*)D^*)K{ID*+X){D^it+) 


0.0006±0.0001 


0.0013±0.0002 


Di*^D(*^X{W*+X){D^Tr+) 


0.00028±0.00006 


0.00021±0.00007 


T-uD(*'**\W*+X){D'^Tr+) 


0.002±0.001 


0.0024±0.0009 


5° ^ i-uD'+{2535){K°D*+){D'^TT+) 


0.553±0.010 


0.620±0.010 


e-i?D*+{2573){K^D*+){D^Tr+) 


0.525±0.010 


0.593±0.009 


D(*W(*')K{£D*+X){D^7r+) 


0.0007±0.0003 


0.0022±0.0003 


Di*^D(*^X{£D*+X){D^Tr+) 


0.00024±0.00008 


0.0012±0.0002 



XII. In the case of the A° baryon, where a fixed sample composition is used, the predicted 
number of events is given as 

Ar(rA+)predicted = Ar(B°)-4v fl + T^l X ^(^° ^ r p,A+X) B(A+ ^ pX-7r+) 

Ju + Jd \ id) 

X£(A°^rP,A+X). (13) 

In this case, the excited charm baryon branching fractions into the A^ state are included in 
the efficiency £(A° — > f^P^A+X) in Eq. (13), as listed in Table Xlll. 

There are four free parameters in the fit for the fragmentation fractions: three relative 
fragmentation parameters /„/ /c^, /s/ {fu + fd), fAt/ ifu + fd), and the normalization parameter 
N{B^), plus three constrained parameters: F, F*, and F**. The values of fg and are fit 
relative to (/„ + fd) to minimize as many biases in the measurement as possible and to high- 
hght the fact that the is reconstructed relative to the and B~ signals. Additionally, 
performing the fit relative to /„ + fa limits any possible inaccuracies in the separation of B^ 
and B~ through the sample composition procedure into the £~D^ and £~D'^ final states. As 
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TABLE XII: Efficiencies in 



'Df Monte Carlo simulation. 



Decays 



e+SVT 



;U+SVT 



£-uD*+{Dtj) 

e.-i>Dt,{2m){7:^D*+){Dti^^) 
i-uD+{2m){D+j) 

ruD*+n^{Dt7) (NR) 
i-i?D+7T^ (NR) 
Dl*^D'^*^X{£DfX) 
Di*)+Di*^~{eD+X) 



T-vD 



{tD+X) 



0.998 • 10-3 (9,980±100) 
1.035±0.014 
0.684±0.011 
0.709±0.011 
0.710±0.011 
0.436±0.008 
0.479±0.008 
0.0023±0.0003 
0.0075±0.0003 
0.034±0.005 



1.201 • 10-3 (12,010±110) 
1.11±0.014 

0.773±0.011 

0.786±0.011 

0.781±0.011 

0.591±0.009 

0.722±0.010 
0.0086±0.0007 
0.0175±0.0005 

0.052±0.001 



B+ D^^^ D^*^- {IDt X) 



0.0055±0.0003 



0.0126±0.0005 



0.0055±0.0003 



0.0109±0.0005 



TABLE XIII: Efficiencies in l~K+ Monte Carlo simulation. 



Decays 




e+SVT 


/x+SVT 






0.629 • 10-3 (6,290±90) 


0.722 • 10-3 (7,220±100) 




-Uh+X{ik+X) 


0.026±0.001 


0.033±0.0007 



mentioned earlier, the sum of fragmentation fractions /„, /d, /s, and /a^, is not constrained 
to unity in the fit, since not all h baryons are necessarily accounted for by reconstructing 
^-A+ states. 

The electron and muon samples are fit separately, since the relative lepton efficiencies 
between the electron and muon modes are not expected to readily cancel. The fit re- 
sults are given in Tabic XIV with statistical errors only indicating good agreement be- 
tween the e+SVT and ^+SVT data sets. As mentioned previously, fq indicates the frag- 
mentation fraction integrated above the momentum threshold of sensitivity in the data, 
fq = fq[prp[B) > 7 GeV/c). Note that about 90% of the i'^A^ combinations in data have 
transverse momenta below ~ 20 GeV/c, but none have pt less than 7 GeV/c. 



VI. SYSTEMATIC UNCERTAINTIES 



The main uncertainties in the measurement of the relative fragmentation fractions come 
from the uncertainties in the branching fractions of the charm mesons, which contribute 
both directly and indirectly to the measurement, and the uncertainty associated with the 



39 



TABLE XIV: Fit results with statistical errors only. 



Fit Parameter e+SVT ^+SVT 

U/fd 1.044 ± 0.028 1.062 ± 0.024 

fs/ifu + fd) 0.162 ±0.008 0.158 ±0.006 

fAjifu + fd) 0.292 ± 0.020 0.275 ± 0.015 

r [ps"i] 0.0157 ±0.0007 0.0154 ± 0.0007 

r* [ps-1] 0.0327 ± 0.0014 0.0331 ± 0.0013 

r** [ps-1] 0.0145 ± 0.0010 0.0146 ± 0.0010 
iV(^O) (10^) 2.02 ±0.07 2.93 ±0.10 



baryon spectrum, which affects the simulation-l:)ased efficiency. The uncertainties in 
the measurement due to the XFT and dE/dx efficiencies are negligible in comparison with 
other systematic uncertainties. The complete list of systematic uncertainties assigned to the 
fragmentation fractions is given in Table XVII. A weighted average between the e±SVT and 
/x+SVT samples is calculated before and after applying a particular systematic variation in 
order to determine the systematic uncertainty for a given quantity. The determination of 
the individual systematic uncertainties is discussed in the following sections. 

A. False Lepton Backgrounds 

The wrong sign lepton-charm combinations represent several possible backgrounds that 
may be present in the right sign signals with a significant contribution to the wrong sign 
combinations expected to arise from false lepton candidates. Another contribution origi- 
nates from real leptons from non-bottom sources, such as electrons from photon conversion 
7 e^e~ or muons from kaon and pion decay- in- flight. These sources are included in 
the discussion of false lepton backgrounds. The wrong sign signals are present in the data 
even after the prompt region is removed by requiring ct*{i~D) > 200 fim, as can be seen in 
Figs. 14 and 15. Additionally, some discrepancy is still observed in the ct*{i~D) comparisons 
between data and simulation, possibly indicating a residual background from false leptons. 
False leptons that originate from a "5"-like hadron {i.e. a relatively long-lived particle) are 
not necessarily represented equally between right sign and wrong sign combinations, as is 
the case with prompt false leptons. Since the false leptons of concern most likely come from 
a real bottom hadron in which a hadronic track has been mis-identified as a lepton, they 
are enhanced in the right sign over wrong sign lepton-charm combinations. 

This systematic uncertainty is studied by utilizing the large false lepton sample available 
from the TTT semileptonic B decays, which has approximately five times more and 
£~D^ events than the ^+SVT trigger sample. Since the statistics are much larger in the 
TTT sample and the average lepton transverse momentum is lower, a larger sample of false 
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FIG. 14: yU+SVT wrong sign invariant mass distributions of (a) Z)°, (b) D*+, (c) D+, (d) Df, and 
(e) A+ after the ct*{fi~D) requirement. 



leptons is available for the study. Lepton candidates with a low probability of being true 
leptons, as measured from a likelihood weighting of lepton identification variables [55, 56], 
are selected from semileptonic TTT events in which one of the charm daughters is matched 
to one SVT track and the lepton is matched to the other SVT trigger track. This sample 
of false leptons is then used to estimate the factor required to scale the residual wrong sign 
"signals" in the £+SVT data to the right sign signals, giving an estimate of the false lepton 
contamination from bottom hadrons in the right sign signals. All selection requirements 
used in this analysis are applied and the numbers of right sign and wrong sign events are 
compared. The scaling obtained for false leptons from long-lived "bottom" -like hadrons is 
We = 2.93 ± 0.47 and = 3.91 ± 0.73, assuming that all of the wrong sign events originate 
from long-lived sources. To obtain the systematic uncertainty, the right sign lepton-charm 
yields are decreased by the wrong sign yields scaled by the appropriate factor for the e+SVT 
and /i-|-SVT datasets. The fit for the fragmentation fractions is then repeated and the 
resulting systematic uncertainties noted in Table XVII. 
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FIG. 15: e+SVT wrong sign invariant mass distributions of (a) D^, (b) D*^, (c) D^, (d) Df, and 
(e) A+ after the ct*{e~D) requirement. 

B. Variation of Selection Requirements 

The selection requirements have been chosen to be similar across the five lepton-charm 
channels, in order to cancel as many systematic uncertainties as possible while still respecting 
the different kinematic features of the decays. To check the dependence of the final result 
on these selection criteria, the signal selection has been varied in such a way that the 
Pt{D) > 5 GeV/c cut is apphed to all channels, while the Xxy{D) and vertex probability 
requirements are the same. The varied cuts used to assign the systematic uncertainty are 
listed in Table XV. 



C. D+ Reflection in the D+ Signal 

In addition to residual wrong sign backgrounds, another source of irreducible non- 
combinatoric background arises from the K^K^n^ reflection into the — »• 

K^Ti^'K^ signal. This effect has been measured from the data, using the simulation to 
scale the expected rates of generic K'K^n'^ decays to the 07r+ decay (see 
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TABLE XV: Alternative signal selection. 
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Section II C 5). By default, the Df reflection is included in the fit to the D'^ signal by con- 
straining the normalization of the reflection within its uncertainty. In order to assign a 
systematic uncertainty on this method, the normalization of the reflection is allowed to vary, 
both by fixing N{Df) to the number measured from data, and also by allowing N{Df) to 
float in the fit to obtain the signal, listed for both scenarios in Table XVI. The larger 
effect is observed when N{Df) is a free fit parameter, while fixing the normalization pro- 
duces a more moderate shift. Since the normalization procedure of N{Df — * </)7r+) relative 
to A^(-D+ — > K~^K~7i^) is, in principle, well-understood from the data and simulation, the 
variations obtained from fixing the normalization are taken as the systematic uncertainty 
associated with this method. 
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D. XFT Efficiencies 



Knowledge of efficiencies that are different for the different particle species is essential for 
the proper determination of the relative efficiencies between lepton-charm channels. One of 
these sets of efficiencies is the XFT trigger efficiencies (described in Section IV Bl), which 
cannot be accurately predicted by the simulation and are not expected to readily cancel in 
the relative efficiencies between the final state charm signals. The systematic uncertainty 
on this efficiency is determined by varying the default XFT efficiencies by 1 u of the fit 
parameters given in Table VIII. To determine the systematic uncertainty of the shift in 
the K- and tt-XFT efficiencies, the tt efficiency, which has the larger uncertainty of the 
two, is shifted up or down by the uncertainties in the XFT parameterizations, while the K 
efficiency is held constant. The systematic uncertainty associated with the proton efficiency 
is assigned by fitting the efficiency with a constant line. The proton parameterization was 
shifted by the full uncertainty on the fit parameters (either all up or all down), which are 
quite large as can be seen in Table VIII. 

E. Single Track Efficiency 

The efficiency to add a single track, needed to adjust the two track topology to the three 
track topology, is measured by reconstructing the K~ti^ti~ti^ channel relative to the 

— » K~TT^ decay mode (sec Section IV B 3). This method assumes that the two additional 
pions in the K^n^n^n^ decay arc uncorrclatcd. Since the two tracks are identified 

and pass through different parts of the detector, this assumption is reasonable. The only way 
the determination of the single track efficiency might be biased by the correlation of the third 
and fourth tracks arises from vertexing effects. To assess the degree of the bias that might 
occur in the vertexing of the K'n'^n'Tr'^ mode due to the correlation between the 

two additional pions, a three track vertex is formed in the Monte Carlo simulation of ii~D^ 
decays with — > K'tt'^tt'tt'^ and the impact parameter of the fourth track with respect to 
the three track vertex is determined. If the impact parameter of the fourth track lies outside 
of 1 (T of the error on the vertex and 1 cr of the error on the impact parameter, it is assumed 
that the fourth track could bias the position of the vertex. The impact parameter of the 
fourth track is found to be outside 1 cr of the three-track vertex (5.3±0.1)% of the time. This 
fraction is assumed to correspond, to good approximation, to the degree of correlation in the 
efficiency of the third and fourth pions. Another source of uncertainty in the determination 
of the single track efficiency is the error on the ratio of branching fractions of the two 
reconstructed decays, B{D^ K-7i+7i-7i+)/B{D^ K-n+) = 2.10 ± 0.03 ± 0.06 [54]. 
The systematic uncertainty from the error on this ratio of branching fractions is also included 
in the systematic uncertainty from the single track efficiency listed in Table XVII. 
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F. Sample Composition Lifetimes 

B meson lifetimes relative to the lifetime are included in the sample composition 
procedure. Consequently, the B lifetimes are needed to determine the predicted number of 
lepton-charm mesons (see Section V). As there are uncertainties on the PDG values [42] of 
the lifetimes used in the fit, which are listed in Table V, the lifetimes and lifetime ratios are 
varied in the process extracting the sample composition within their PDG uncertainties. The 
central value of the hfetime ratio t{B~)/t{B'^) — 1.086 ±0.017 has changed several times in 
several years and different values are used in the sundry measurements of fu/ fa- Although 
the hfetime ratio in the PDG is shghtly higher than that used in other measurements of /„/ /d, 
the uncertainty on the PDG value covers the central value of the other possible hfetime ratios. 
No A° lifetime is used in the fit for the baryon sample composition, although f^^/ (/„ + f^) 
varies shghtly when the ratio t{B~)/t{B^) is varied within the PDG uncertainty, because 
the A^ fragmentation fraction is measured relative to the B^ and B~ modes. 

G. Monte Carlo Simulation Statistics 

Since a finite number of Monte Carlo simulation events are generated for each exclusive 
decay to be used in the process to extract the sample composition (see Section III,) the 
statistics of the generated simulation is checked to see whether the statistical uncertainties 
on the yields, which are used to determine the efficiencies, contribute a significant uncertainty 
to the measurement. The simulation yields in each decay arc shifted by ± 1 o" around their 
central values and the efficiencies arc rc-dctcrmincd accordingly. To assign the systematic 
uncertainty, half of the yields are randomly shifted up, while the other half are shifted down. 
In all cases, the shift in all three relative fragmentation fractions is small compared to the 
other uncertainties, as can be seen in Table XVII. 

H. Bottom Hadron Lifetimes 

Knowledge of the bottom hadron lifetimes is also needed for the generation of the various 
B simulation samples. While the B^ and B~ lifetimes are well-measured, there are large 
uncertainties on the S° and A° hfetimes [42]. To assign a systematic error due to the 
uncertainty in the knowledge of the B and A^ lifetimes, the simulation is re-generated with 
the B^g and A° lifetimes shifted by one sigma uncertainty on their PDG values: r(5°) = 
(438 ± 17) /im and r(AO) = (368 ± 24) /xm. The shift in fxjifu + fd) is one of the larger 
uncertainties in Table XVII, but it is still small compared to the uncertainties due to the 
imprecise knowledge of the baryon branching fractions. 
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I. pt Spectra 

The bottom hadron pt spectra are one of the biggest uncertainties on the knowledge 
of the relative efficiencies. Consequently, the systematic uncertainties arising from the px 
spectra are estimated conservatively, as no definitive measurements for the meson and 
A° baryon are available. The systematic uncertainty assigned to the pt spectrum for £~A+ 
is taken from a ±2 a variation of the tuned semileptonic A° px spectrum. The spectrum 
is varied by ±2 cr in order to provide a conservative error, since the pr spectrum measured 
from the semileptonic A° decay is incomplete. Although the l~D'^ MC generated with the 
Pt spectrum obtained from the inclusive J/%1) cross-section measurement [20] agrees well with 
the data, there is the possibihty that the S° meson pt spectrum is different from the and 
B~ spectra. This possibility is accounted for by measuring the ratio /«/ (/„ + fd) with the 
default generator input spectrum, while the S° and B~ decays are generated with the px 
spectrum inferred from Ref. [20]. Since no significant discrepancy is observed between the 
£~D^ data and the sinnilation using the pt spectrum from Ref. [20], this is a conservative 
assessment of the systematic error due to the uncertainty on the B^ momentum spectrum. 

J. Specific Ionization Efficiency 

Accurate knowledge of the requirement on the dE/dx based likelihood ratio CTZ{p) on 
the proton in the A+ pK~Ti^ decays is important for an accurate determination of the 
^^A+ efficiency relative to the semileptonic bottom hadron decay efficiencies. In order to 
assign an uncertainty to the knowledge of the dE/dx efficiency, the measurement of the 
fragmentation fractions is performed without any dE/dx cut applied to either the data or 
simulation. The fA^/{fu + fd) fit result is stable and the difference with the default fit is 
treated as a systematic uncertainty. Removing the dE/dx cut does not produce a significant 
change in either or fs/{fu + fd)- 

K. Polarization 

The polarization of the A^ baryon in hadronic collisions is not known. By default, the 
A° baryon is unpolarized in the simulation used in this measurement. In order to assign a 
systematic uncertainty to the possible polarizations of the A°, the extreme cases of the A° 
being fully polarized are tested to bound the effect. A systematic uncertainty is assigned 
when the A^ is produced either with entirely positive helicity or entirely negative helicity. 
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L. Bottom Hadron Branching Fractions 

Systematic uncertainties due to the knowledge of the bottom hadron branching fractions 
arise in two places in the sample composition procedure. First, the indirect semileptonic de- 
cays {e.g. B ^ DDX, D £^X) contributing to the lepton-charm signals, many of which 
are poorly determined experimentally, and second the uncertainty in the PDG semileptonic 
A° branching fraction, i3(A° — > £~z/^A+X) = (9.2 ± 2.1)% [42]. Since many of the measured 
indirect semileptonic decays are poorly determined, the branching fractions predicted from 
symmetry principles for these decay modes are used in the sample composition process to 
determine the systematic shift in the fragmentation fractions. The contributions of indirect 
semileptonic decays to the i~D'^, £~D^, and l~D*^ signatures are small (see Tables IX- 
XI), but the rate for B^/B' D^*^'^ D^*\Wf X) decays contributing to the i'Df final 
state (see Table XII) is an order of magnitude larger due to the more copious fragmentation 
of b quarks into and B~ mesons versus B^ mesons. To give a sense of the maximal possi- 
ble effect on fs/ {fu + fd) if no contributions from B^ and B^ mesons to the i^Df yield were 
accounted for, the fragmentation fraction fs/ {fu + fd) would increase by about 10% from 
~ 0.160 to ~ 0.176. This estimate is presented for general interest, though it not used as a 
systematic uncertainty on the measurement, as it is known to be an incorrect assumption. 

To determine the systematic uncertainty associated with the inclusive semileptonic 
A° branching fraction, the PDG value is varied within its uncertainties. This is one of the 
largest systematic uncertainties associated with the measurement of the A^ fragmentation 
fraction. 

M. Charm Branching Fractions 

Another source of systematic uncertainty due to the branching fractions used in the 
sample composition procedure arises from the often poor knowledge of the ground state 
charm branching fractions, which are taken from the PDG and listed in Table VII. To 
determine the uncertainty in the fragmentation fractions, the central values of the ground 
state charm branching fractions included in the sample composition are varied, one by one, 
within ±1(1 of the PDG uncertainty. The largest shift in fu/fd comes from B(D'^ — > 
K~7r^7r~^), while the single largest uncertainty in /«/(/« + fd) is due to the large error 
on B{D^ 07r"^). A poor knowledge of -B(A+ pK'-K^) contributes the largest single 
systematic uncertainty to /a(,/ {fu + fd)- 

In addition to the poorly measured ground state charm branching fractions, many of the 
excited charm decays also have large uncertainties. To assess a systematic uncertainty for 
the limited knowledge of the excited charm decays, the excited charm branching fractions are 
varied by shifting half of the D** branching fractions randomly up by 30%, while the other 
half are shifted down. When quoting the final result on the fragmentation fractions, a sep- 
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arate systematic uncertainty is quoted due to uncertainties on external branching fractions 
as indicated in Table XVII. 



N. Sample Composition 

A systematic uncertainty is assigned for the uncertain knowledge of the sample com- 
position (see Section 111). The sample composition is considered without any of the 
non-resonant baryon modes included, while the total semileptonic branching fraction in 
both cases is required to be B{A^ i^UiA'^X) = 9.2%. A systematic effect for a potential 
mis-modeling of the decay is also considered and found to be negligible. This uncertainty is 
determined by evaluating the width difference of the m(£~A+) distribution in both the data 
and the Monte Carlo simulation. The RMS of the data distribution is 451 MeV/c^, while 
the RMS of the MC is 455 MeV/c^ The ratio of excited to ground state A^^ i~i^iA+ de- 
cays is changed in the simulation such that the RMS of the simulated m(£~A+) distribution 
decreases by 4 MeV/c^, producing a 0.017 shift in fKb/ifu + fd)- Both uncertainties result 
in a total systematic uncertainty of 0.047 on /aj/(/u + fd)- 



O. Total Systematic Uncertainty 



The total systematic uncertainties due to the knowledge of the relative efficiencies, ob- 
tained by adding the individual systematic uncertainties related to the determination of 

the efficiencies used in the parameterization of the sample composition in quadrature, are 
{laol} for /„//,, for /./(/„ + /,), and for .f.^J{k + h). When uncertain- 

ties arising from branching fractions are included, the uncertainties increase to {^0 074} for 
fu/fd, {ISiSi} for fs/Uu + /d), and for UJ^fu + /d), as given in Table XVII. 



VII. FINAL RESULTS AND DISCUSSION 



The weighted average of the fragmentation fractions between the electron plus displaced 
track and muon plus displaced track samples yields: 

^ = 1.054 ± 0.018 (stat) 1°-°^^ (sys) ± 0.058(B), 

Id 

/•S I A AAr +0.011 +0.057 



= 0.160 ± 0.005 (Stat) t'^o^ (sys) t'^^ (B), 



fu + fd 

. = 0.281 ±0.012 (Stat) 1°:°^^ (sys) 1°:J^?(S). 



fu + fd 

Since this analysis potentially ignores the production of b baryons that might not decay into 
the A° final state, no constraint is applied requiring the fragmentation fractions /„, fa, fs, 
and to sum to unity. The correlation matrix for the fit is shown in Table XVIII. 
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TABLE XVII: Compilation of systematic uncertainties assigned. 
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This result is in agreement with the world average of the fragmentation fraction of 
B~ relative to which is expected to be equal to unity [42]. The result on the rela- 
tive fragmentation fraction fg/ (/„ + fa) presented in this paper agrees with the LEP average 
fs/ifu + fd) = 0.135 ± 0.011 [42] within one standard deviation. Separating B{Df — > ^tt"*") 
from the result of fg/ {fu + fd) , for comparison with the world average, gives: 

X B{D+ 07r+) = (5.76 ± 0.18 (stat) toil (sys) ) x 10"^ 

Ju + Jd 

There is no significant indication of a higher rate of 6-quark fragmentation to mesons 
at the Tevatron which would contribute to the anomalous Run I values of % [11, 12, 19]. 
The uncertainty on fg will significantly decrease with an improved measurement of B{Df —>■ 
077"*"), which is in preparation by the CLEO-c experiment [57]. 

Separating the poorly known B(A+ — > pK'i:^) and -B(A° — > .£~P^A+) from the results of 
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TABLE XVIII: Correlation matrix of fit parameters. 
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-0.575 


r 




1.0 


0.657 


-0.122 
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fhjifu + fd) yields: 

Ju + Jd 


>pK-T:+) = (14.1 ± 0.6 (stat) 


(sys) ) X 10- 


or 



X B{hfi rPeAt) B(A+ ^ pK-TT+) = (12.9 ± 0.6 (stat) ± 3.4 (sys) ) x lO'l 

fu + fd 

This quantity can be compared more naturally with the LEP results, which quote f\^ x 
i3(A^ £"z/^A+) X i37^(A+ — pK~7r~^) [58, 59]. When all branching fractions with large 
uncertainties are factored out, /a(, is ~ 2.3 ex higher than the LEP results, assuming that 
fu = fd = 39.7% at LEP. In addition, this measurement of fKtKfu + fd) is approximately 
twice as large as the world average of /aj/(/« + fd) = 0.125 ± 0.020 [42], which is dominated 
by the LEP results. 

A. Discussion of Results 

The difference between the /a^ result presented in this paper and the LEP results may 
be explained, at least in part, by the different environment of hadro-production of bottom 
hadrons in pp collisions. In addition to this effect, the transverse momentum of the bottom 
hadrons is significantly lower for the data used in this measurement, {prip) ) ~ 15 GeV/c, 
than the hh data collected at the Z pole used in the LEP measurements, (prib) ) ~ 45 GeV/c. 
To study a potential momentum dependence of /aj,, the behavior of the fragmentation frac- 
tions in bins of the lepton-charm p^ is investigated. Note that the fragmentation fractions 
can depend on momentum, and the fractions reported here are for momenta integrated above 
the effective p™" which is chosen to be 7 GeV/c in this analysis. 

For this study the electron and muon datasets of the lepton-charm candidates are divided 
into three momentum ranges with similar statistics in each bin. The chosen momentum bins 
are less than 11 GeV/c, from 11 to 14 GeV/c, and greater than 14 GeV/c. The lepton-charm 
yields and corresponding efficiencies are redetermined in each momentum interval and the fit 
for the fragmentation fractions is repeated. The weighted average of /a^/ (/„ -l- fd) obtained 
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FIG. 16: Determination of /a^/ (/u + fd) for three momentum ranges (three points with error bars) 
overlaid on Monte Carlo simulation scaling of /a^/(/u + fd) as a function of p'r(^~A+) (triangles). 
The dashed line is a fit to the Monte Carlo simulation. 



from the e+SVT and /i+SVT data in the three pt ranges is shown as three points with 
error bars in Fig. 16. The uncertainties on these points include the systematic uncertainties 
on the efficiencies, but do not reflect the uncertainties from branching fractions. The data 
points are consistent with a decrease in the ratio fAi/{fu + fd) with increasing bottom hadron 
momentum. However, in the near future, larger CDF datasets of lepton-charm events will 
provide increased statistics for a more adequate extrapolation of this suggested momentum 
dependence of /aj, as compared to S-hadron momenta at LEP. 

To obtain a better extrapolation of the indicated momentum dependence of 6-quark frag- 
mentation into baryons, the Monte Carlo simulation tuned on the data is used to estimate 
such a momentum dependence. The inclusive — £~i/£A+X simulation generated with 
the tuned semileptonic A° spectrum is compared with the inclusive I'vpD^X 
simulation generated with the pt spectrum inferred from the inclusive J/ip cross-section 
measurement [20]. Assuming that the inclusive Monte Carlo samples provide a good de- 
scription of the ^~A^ and i~D~^ data, as demonstrated in Section IV A 1, the ratio of the 
momentum dependence of both datasets provides an estimate of the shape of the ratio of 
dA^/dpr(A°) to dN/dpT{B^)- The ratio of both distributions is therefore proportional to 
fki/fd and thus to fA^/{fu + fd), assuming = fd. To obtain an absolute normalization, 
the ratio of the two distributions is fixed at the mean pT{i~A^) of the present measure- 
ment, (pt(^~A+) ) ^ 14.1 GeV/c, to the central value of fh^/{fu + fd) = 0.281 as obtained 
in this analysis. The result of this MC study is shown as triangles with error bars in 
Fig. 16. Fitting a straight line to these points agrees well with the three data points of 
fAb/ifu + fd) obtained in three momentum bins as described above. Extrapolating the line 
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to Pt(^"AJ) ~ 35 (40) GeV/c yields a value for /a,/(/« + fd) of about 0.128 (0.092). This 
is close to the world average of fh^l{fu + fa) — 0.125 ± 0.020 [42], which is dominated by 
the LEP results. This study indicates a possible momentum dependence of the 6-baryon 
fragmentation that would explain the difference between the result presented in this 
paper and the LEP measurements. 

Finally, knowledge of /a^ will improve with better measurements of the — >• 
pK~7r~^ branching fraction and the semileptonic — * £^i?eA^X branching fractions, in 
addition to better measurements of the A^ semileptonic sample composition, particularly 
the measurement of the A^ — > £~i/£Ac(2593)+ and A^ — > £^z/f A^(2625)"'" branching fractions. 
Additionally, a definitive measurement of the pt spectrum of baryons in pp collisions com- 
pared to the momentum spectrum of B'^ mesons measured with fully reconstructed A^ and 
S° decay modes will shed light on expected differences in the momentum spectra and sig- 
nificantly reduce the systematic uncertainty of a measurement of the 6-quark fragmentation 
fractions in the future. 
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